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WETTING DYNAMICS OF EVAPORATING DROPS ON VARIOUS SURFACES 


by F. F. Simon and Y. Y. Hsu 
Lewis Research Center 

National Aeronautice and Space Administration 
Cleveland, Ohio 

SUMMARY 

Shape histories of water drops evaporating at room temperature on 
copper, Incite, stainless steel and teflon were recorded. The drops were 
observed to go through three stages. In the first stage, the interface radius 
remained constant while the contact angle reduced with decreasing volume, 
until the receding angle is reached. In the second stage, the radius de- 
creased almost linearly with time while the contact angle remained a little 
less than the receding angle. In the third stage, both radius and contact 
angle reduced rapidly. The relative leng;th of each stage is determined by 
the type of surface. For a highly wettable surface, Stage I dominated, for 
a very nonwettable surface, stage II dominated. An analysis is presented 
enabling the prediction of shape history of a drop provided the evaporation 
rate and flow resistance over the surface are given. 

INTRODUCTION 

The wetting dynamics of an evaporating liquid film is important in de- 
termining the efficiency of heat transfer in two-phase flow and the stability 
of a liquid film in surface coatings. In the annular flow region of two phase 
flow there is contact of a liquid film with a heated surface. As this film 
evaporates fluid into the core, a nonwetted surface condition may develop. 
The growth of nonwetted areas is a limiting condition to good heat transfer 
in the annular flow regime. Since the initiation and growth of a nonwetted 
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condition will be determined to a great extent by the wetting nature of the 
evaporating liquid film on the solid surface, a study of wetting characteris- 
tics during evaporation is usefuh Another area in which wetting is important 
during an evaporation process is the stability of liquid surface coatings. 
Surface coatings such as paints may have a tendency to ’’crawl” and leave 
a portion of the surface uncovered. While this situation is initiated by such 
things as poor surface preparation, evaporation may further increase the 
tendency of the coating to ”crawl”. Therefore, the dynamic wetting of an 
evaporating liquid film as a function of different surface condition is of 
interest to both situations above. 

An index of wetting characteristics is the contact angle. In the case of 
a static interface, the contact angle is related to .the energies at the triple 
interface line of gas, solid and liquid (Fig, 1(a)). This energy relationship 
is stated by the following equation of Young and Dupre. * 

y LV cos 6 e - y gv + y gL = 0 (1) 

In a dynamic situation where the interface is moving, such as when the 
liquid is being removed by evaporation, equation (1) is not expected to hold. 
Under dynamic conditions, the contact angle has to be increased beyond a 
critical value before the interface will move from a wetted to a nonwetted 
area. Conversely, an interface would not recede until the contact angle is 
reduced below a certain angle. These two critical angles are the well-known 
equilibrium advancing and receding angles. The Young- Dupre equation was 

(See Nomenclature for the definition of symbols. ) 
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modified in Ref. 1 to account for the hysteresis of contact angles as follows* 
y L y(cos 6 r ) + Cy SL - Cy gv - F r = 0 for receding angle (2a) 

y LV cos # a + Cyg L - Cy sy + F & = 0 for advacing angle (2b) 

The terms F and F account for the additional surface force per unit length that 

i a 

is required for the interface to begin to move (Fig. 1(b)). The energy in- 
volved is called "Contortional Energy" by Good (Ref. 2). The coefficient 
C is to account for the increase of area due to the presence of roughness. 

The equilibrium advancing and receding angles can be measured by 
tilting a plate to cause the drop to deform until it slides. Another method 
of obtaining an advancing contact angle is by the addition of fluid to a drop. 
Conversely a receding angle can be achieved by removal of fluid from a 
drop. One of the ways of removing liquid is to measure the angle change 
during evaporation, provided the evaporation rate is extremely slow (Refs. 3 
and 4). In Ref. 4, the receding angle was measured by allowing a drop to 
evaporate and noting the condition at which the drop radius at the solid- 
liquid interface decreased for constant contact angle. Of special interest 
to this study are the three stages noted in Ref. 4 which occurred during drop 
evaporation. At the initiation of drop evaporation the contact angle decreased 
with time, while the solid-liquid interface radius remained constant (stage I). 

Stage II became evident when the contact angle remained approximately con- 
stant with time (receding contact angle) and the solid liquid interface radius 
decreased. Stage III occurred at the end of the drop life when both the con- 
tact angle and the solid-liquid interface radius decreased with time. 

It is the intent of this paper to indicate the variation in interfacial be- 
havior that can be expected from a receding triple- phase line under condi- 
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ti oiis of evaporation. The contact angles of evaporating water drops on four 
surfaces with a wide range of wettability were measured. The measured 
histories of contact angle and the liquid-solid interface radius were compared 

with an analysis. 

EXPERIMENTAL METHOD AND PROCEDURE 

The basic experimental approach was to record the contact angle and 
size of an evaporating water drop as a function of time on surfaces of teflon, 
Incite, copper and stainless steel. Drops were evaporated naturally under 
ambient conditions. Figure 2 shows the position of the drop with respect to 
movie camera and to a microscope with a goniometer. The camera was run 
at a speed of 1 frame every 11 seconds and was synchronized with two elec- 
tronic light sources. A millimeter scale was placed immediately in back 
of the drop to provide a length reference. The microscope with the goniom- 
eter attachment was used to measure contact angles of the evaporating drops 
periodically. Although the views of camera and goniometer are perpendicu- 
lar to each other, the near symmetry of a drop minimized the possible dif- 
ference in contact angle between the side and front views. 

Prior to the runs using copper or stainless substrates, the surfaces 
were polished. A metallographic polishing procedure was performed which 
went from a rough grind to a final step of polishing in a vibratory polishing 
machine with levitated alumina grit. All surfaces were ultras onically 
cleaned with mild detergent, rinsed with distilled water and dried with 
filtered air. The procedures stated above for sample-preparation were 
followed each time to ensure reproducibility of results from run to run. 
Reproducibility was important in making independent measurements of re- 
ceding ingle which serve to characterize the surface. 
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The 16 mm movie films of the evaporating drops were placed in an 
analyzer that permitted measurement of the projected image. In this way 
it was possible to obtain the liquid-solid interface radius as a function of 
time. By using the radius and contact angle information, calculations were 
made of the volume and evaporation rate. 

For the measurements of equilibrium receding angles, the sliding-drop 
method was used. The deformation and motion of a drop on a slowly tilted 
surface was recorded by motion pictures from which the receding angles 
were measured. 

RESULTS 
Experimental Data 

Typical curves of liquid-solid interface radius and contact angle against 
time during drop evaporation are shown in Figs. 3-6. These figures are 
arranged in the order of increasing receding contact angles; namely, water- 
copper, water-lucite, water- stainless steel, and water-teflon. With ex- 
ceptions to be discussed later, these figures have features common to those 
discussed in Ref. 2. An idealized model forthe general behavior is depicted 
in Fig. 7. First, in the early stages of evaporation there is a constant radius 
period in conjunction with a decreasing contact angle (Stage I). Then, one 
may observe in some cases a second stage (Stage II) in which the contact 
angle is nearly constant but the radius is decreasing almost linearly with 
time. Then there is the last stage (Stage III) in which both radius and con- 
tact angle decrease rather rapidly with time. The lengths of the various 
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stages depend upon the rate of evaporation and the size of the receding angle. 
All the surfaces studied except copper exhibited the three stages indicated 
above. A water drop on copper exhibited Stage I for the majority of the drop 
life. Representative photographs from the movie film of the teflon and copper 
surfaces are shown in Figs. 8 and 9. These pictures demonstrate the fact 
that a teflon surface permits a long Stage- II event while an evaporating drop 
on copper is primarily a Stage- 1 phenomenon. 

Analysis 

In order to describe the variation of drop size with time during evapora- 
tion, a model is proposed as follows: 

(1) The shape of a sessile drop is really not a truncated sphere and can- 
not be stated in explicit form. However, for simplicity a truncated sphere 
geometry is assumed to represent the drop shape. (See Fig. 1. ) 

V = - R 3 (l - cos e) 2 {2 + cos e) = - r 3 l 1 -- - CQS 9 ) i 2 + C . 9 S ^ ( 3 ) 

3 3 3 


with 

R = — ~ — (4) 

sin 6 

It is found that by comparing the true volume of a drop with that for a trun- 
cated sphere, an error of about 5 percent results. 

(2) The drop volume is reduced due to evaporation by the rate 


— = -FS 
dt 


(5) 


with the drop surface S defined as 


27rr 2 

• 2 „ 
sm 6 


(1 - cos 6) 


( 6 ) 
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When the drop volume is reduced, the volume change has to be accomplished 
by either reduction of angle or radius, or both: 


dv _ av 

d 6 + dV 

dr 

(7) 

dt 3 6 

r dt Sr 

e dt 

_rs = — 

1 d0 + SV I dr 

( 8 ) 

36 

r dt Sr 

fldt 



Using equation (8) and the partial derivatives from equation (3), 
the following relation between r and 6 as function of time t can be ob- 
tained . 


- 2r= (l- cos 6) f (2 + cos e) sin + r de' 

siu 2 6 L dt dtj 

The solution of equation (9) depends on the relation between r and 
can be described by the stages of evaporation: 

(a) Stage I: 


(9) 

6 which 


when 


and 


dr 

dt 


0 


e > e 


r J 


( 10 ) 


r = const = r (10a) 

o 

Thus the solution is simple 

(b) Stage II and III: The relation between r and 6 is supplied by a 
force balance to be discussed as follows. 
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Recall in equation (2) that for the equilibrium receding angle, a surface 
force F r was required for an interface to begin to move. When the drop 
is evaporating and its radius is reducing at a finite rate, there is an actual 
flow of liquid which is more than the infinitesinal rate implied in the re- 
ceding angle equations. An additional force is required to effect the flow. 
Although it is theoretically possible to determine the shear stress required 
to cause this flow by considering the basic momentum and continuity equa- 
tions, if is assumed that the additional force required is proportional to the 
time rate of change of the liquid-solid interface radius (Fig. 10). Essentially, 
this model assumes a linear velocity profile or an analogy to the Couette 
type of flow. Figure 10 can be compared with Fig. 1(a) and 1(b). 

F.= -f— (11) 

1 dt 


or, the force-balance equation is modified to read: 


^LV cos 


+ Cy LS = Cy, 


SV +F r 


dr 

dt 


( 12 ) 


Note that when dr/dt - 0, equation (12) is reduced to equation (2) which 
is the case when the interface just begins to move. Substruction of equa- 
tion (2) from the above equation results in: 

y T -rr(cos 6 - cos 6 ) = -f — (13) 

ijV dt 

Differentiation of equation (13) yields: 


y LV sin 


d6 

dt 



(14) 



thus 
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dfl_ f d 2 r 
dt y LV sin 6 ^-2 

with the boundary conditions: 


r = r 


o’ 


dr 

dt 


0 


( 15 } 


at 


0= e r 

Equation (15) is the relation between r and 9 needed for the completion of 
the governing equations 0 

Elimination of d0/dt between equation (9) and (15) yields: 


2 

-2T sin^0 = — — (1 - cos 6) r ^ + (1 - cos 6 )( 2 + cos 0)sin* 


y 


LV 


dt 


dr 

dt 


The initial condition for Stage-II is 


t = t 


II, o’ 


r = r 


o’ 


II, 0 


= 0 , 


(17) 


Thus a complete solution for r and 6 in Stage-II requires the simultaneous 
solution of equations (16), (3), and (5)„ However, since equation (16) is a 
function of r and 6 , it appears that equation (16) and (17) can be used to 
solve for r(t) if 0(t) is known „ Thus, if a proper 6(t) curve can be assumed 
an approximation for r(t) can be obtained which then can be used as an inp 


to an iterative scheme for the complete solution. 



10 


Approximate Solutions 

The inputs needed for solution of equation (16) are: 

(1) Evaporation Coefficient T 

(2) Resistance Coefficient f 

(3) Time history of contact angle 0(t) 

(4) Initial slope of radius r’ 

II,o 

Of the four items, V and f are parameters. Theoretically, they should 
be obtainable through basic consideration of energy and momentum trans- 
port. For lack of simple solutions, we can use empirical values obtained 
from experimental results of r(t) and V(t), with the understanding that 
eventually they should be obtained a priori. It should be noted that although 
empirical values of f are fairly constant over the majority of Stage II, the 
empirical evaporation coefficient r does vary by a factor of 2 with drop 
size. Therefore there is the option of using a mean value of T, (i. e. , r mean ), 
or using T as function of time T(t). 

The third item of the input, 0(t), is the functional relation of contact 
angle with time. It can be either the result of a complete solution, as 
mentioned in the last section, or an approximate form can be used as input 
for iteration. There are two possible choices. First is to use the experi- 
mentally determined history of 6, i. e. , 6 (t). The second possibility to 

assume 6 = 6^ since this is approximately the case for the main part of 
Stage II, 
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Item four is the initial value of r Q at the beginning of Stage II, Theo- 
retically, rl should be zero but experimentally, it is finite. 

Comparison of Experimental Results with Approximate Solution s 
For the comparison of experimental results with analysis, we can ignore 
Stage I since it is simply a constant radius region. The important test is 
after the beginning of Stage II, 

In order to check the analysis and its approximations for Stages II and 
III, machine computations were made using equation (16). The analytical 
results of Fig, 4(a) were made by using experimental values of 0(t)', r^, f 
and T(t) as input to obtain radius against time curves. The agreement of 
the computed values with the experimental data shows that the analytical 


model is realistic. Additional simplifications can be introduced to make 
the model more practical. To show the effect of assuming r^ = 0 at the 

^ jr 

beginning of Stage II, two curves using r’ Q = 0 and r^ = -6.93x10 cm/sec 
(experimental value) were obtained by use of equation (16) with the experi- 
mental value to 6(t) as input. The result is shown in Fig. 6(a). As the re- 
sults show, the assumption of r^ = 0 will introduce an error of about 
5 percent on the average for the life time of the drop. 

To test the assumption of using 6=6 and also r - T three 
curves were generated for the following conditions: 


6=6 (t) 

exp v ' 

r = r exp<« 

r’ = 0 
0 

6=6 (t) 

6Xp v ’ 

r= r „„ 

mean 

r ? = 0 
0 

11 

r = r 

mean 

r’ = 0 
0 


The comparisons are shown in Fig. 5. The comparisons show that the first 
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order approximation of assuming 6 = 6 ^ and r = r mean , = 0 does not 
cause great deviation from the experimental data. 

ADDITIONAL OBSERVATIONS AND DISCUSSION 

In Stage II, the contact angle departs slightly from its equilibrium re- 
ceding value prior to the movement of the triple interface line. This de- 
parture is not always small as can be seen in the curve for teflon (Fig. 6) 
at a time of 1000 seconds. The teflon curve also dramatizes the recovery 
back to the receding contact angle. Apparently, the liquid drop is stressed 
prior to the movement of the triple interface line. In Stage III there is no 
recovery to a receding contact angle. A possible cause for this is that the 
ratio of surface area to volume has become large, causing a more rapid 
decrease in volume than occurred in pahse II. This rapid volume decrease 
does not permit sufficient time for recovery of the drop contact angle. 

The equilibrium receding angle, measured by the sliding-drop method 
are indicated on Fig. 3-6 for comparison with the evaporating drops. With 
the exception of stainless steel, the sliding-drop receding contact angles are 
consistent with the evaporating drop results. For the case of lucite (Fig. 4(b)) 
and teflon (Fig. 6(b)) there is little difference between the receding angle of 
the sliding drop and the Stage-II contact angle. This small difference is to 
be expected because of the very low but finite evaporation rate. Therefore, 
the evaporating drop is a satisfactory method for measuring receding con- 
tact angle provided the evaporation rate is sufficiently low. 
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In our experiments a water drop on a copper surface did not go through 
a Stage II performance. The drop history during Stage II is dependent on 
the receding contact angle. This suggests that the receding contact angle 
for our copper surface was so low that by the time this angle was reached 
the Stage III effects dominated the drop life. The receding angle from sliding- 
drop method verifies this explanation (Fig. 3(b)). 

CONCLUSIONS 

From the experimental results and analysis, the following conclusions 
can be reached. 

1. The change of shape of an evaporating drop is dependent upon the 
type and condition of the solid surface and on the evaporation rate with the 
receding angle as a controlling factor. In the extreme cases, the drop 
radius remains constant while the contact angle decreases continuously for 
a highly wettable surface, while for a highly nonwetting surface, the contact 
angle remains close to the equilibrium receding angle and the radius shrinks. 

2. The shape history of a drop as decribed in Ref. 2 for limited cases 
are found to be applicable to a wide range of surface conditions . It can be 
qualitatively described as composed of three stages: 

Stage I contact angle 6 S: equilibrium receding angle $ r , r = r Q 

Stage II 6 < 6 r , r < r Q 

Stage III 6 and r rapidly approach zero 

3. An analysis has been conducted to predict the shape and size history, 

A first approximation can be obtained by use of a few simplified assumptions. 
Comparison of experimental data with the computations show that these ap- 
proximations permit close agreement with the experimental results. 
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4, It should be cautioned that when the method of evaporating drop is 
used to measure the receding angle, the evaporation rate should be extremely 
slow. Otherwise, an angle smaller than the true receding angle (Stage n) 

will be obtained. 
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APPENDIX - SYMBOLS 
roughness coefficient 

surface force per unit length, or contortional energy , needed 
to start the advancing of an interface 
surface force per unit length resisting the motion of interface 
surface force per unit length, or controtional energy, needed 
to start the receding of an interface 
coefficient associated with 
radius of curvature 
contact radius 
initial contact radius 
= dr/dt 

liquid- vapor surface area of the drop 
time 

time at beginning of Stage II 

volume of the drop 

evaporation coefficient 

mean value of evaporation coefficient 

interfacial tension between liquid and vapor 

interfacial tension between solid and liquid 

interfacial tension between solid and vapor 

contact angle 

advancing contact angle 

static contact angle with interfacial tensions in equilibrium 
receding contact angle 
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(A) GEOMETRICAL REPRESENTATION OF 
A DROP WITH TRIPLE INTERFACE 
W VECTORS. 


VAPOR 



(B) FORCE BALANCE FOR AN EQUILIBRIUM 
RECEDING ANGLE. 

Figure 1 



Figure 2. - Apparatus configuration. 
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Figure 5. - Shape history of a drop, radius and contact angle. 
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(B) 

F igure 6. - Shape history of a drop, radius and contact angl 




Figure 8. - Evaporation of a water drop on a teflon surface. 



t ■ 4050 SEC 

Figure 9. - Evaporation of a water drop on a copper surface. 
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Figure 10. - Force balance on a moving inter 
face (eq. (9)). 



